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ABSTRACT 

Tidal disruption events occur when a star passes too close to a massive black hole and it is totally ripped apart by tidal forces. It may 
also happen that the star is not close enough to the black hole to be totally disrupted and a less dramatic event might happen. If the 
stellar orbit is bound and highly eccentric, just like some stars in the centre of our own Galaxy, repeated flares should occur. When 
the star approaches the black hole tidal radius at periastron, matter might be stripped resulting in lower intensity outbursts recurring 
once every orbital period. We report on Swift observations of a recent bright flare from the galaxy IC 3599 hosting a middle-weight 
black hole, where a possible tidal disruption event was observed in the early 1990s. By light curve modelling and spectral fltting we 
can consistently account for the events as the non-disruptive tidal stripping of a star into a highly eccentric orbit. The recurrence time 
is 9.5 yr. IC 3599 is also known to host a low-luminosity active galactic nucleus. Tidal stripping from this star over several orbital 
passages might be able to spoon-feed also this activity. 
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1. Introduction 

Supermassive black holes reside at the centre of most galaxies. 
When accreting matter, these black holes emit a huge amount 
of energy, becoming active galactic nuclei (AGNs). The energy 
and radiation produced by matter accretion onto these black 
holes play an important role in determining their masses and 
spins as well as the properties of the host galaxy bulges. AGNs 
usually have a relatively short duty cycle and are in a low- 
luminosity state for most of the time (Ho 2008). AGNs may be¬ 
come bright when major inflows of gas feed the compact object 
(Hopkins et al. 2006). Apart from intense episodes of accretion, 
the black hole is powered at a floor minimum level mostly by 
a hot gas from colliding stellar winds in the central region. A 
star in the nuclear cluster can be scattered close to the black hole 
via dynamical encounters and put into an orbit too close to the 
black hole. When the black hole tidal force overcomes the star 
self gravity (at the so called tidal radius) the star starts to be 
ripped apart. About half of the star mass remains bound to the 
black hole, forming an accretion disc and powering a luminous, 
long lasting (months to years) accretion flare (e.g. Rees 1998; 
Phinney 1989). Despite their elusiveness, a number of events 
have been reported in the literature, resulting in the complete dis¬ 
ruption of the star (e.g. Renzini et al. 1995; Bloom et al. 2011; 
Gezari et al. 2012). 

Tidal disruption events (TDEs) provide an unique probe to 
reveal the presence of otherwise quiescent black holes, allowing 
us to study the accretion process also in a different regime than 
that of AGNs (Komossa 2012; Gezari 2012). TDEs were first 
predicted and later on observed in the soft X-ray to ultraviolet 
bands where the peak of the accretion disc emission lies (Strubbe 
& Quataert 2009). Two relativistic TDEs were discovered based 
on high energy triggers, thanks to the launching of relativistic 
jets along the line of sight (Bloom et al. 2011; Burrows et al. 


2011; Cenko et al. 2012). Only a small number of TDEs how¬ 
ever have been observed to date in X-ray, UV and optical sur¬ 
veys, mainly because of their low rate of occurrence and sparse 
observations (Donley et al. 2002). 

TDEs are ultimate events. A star entering the tidal radius rt ~ 
R* (where Mbh is the black hole mass and M* and 

R, are the mass and radius of the star) is completely disrupted 
by tidal forces. The fate of a star orbiting its central black hole 
is then defined by its pericenter radius, r^: if j3 = r^|r^ i 1 we 
have a TDE. What happens instead if fi f 1 ? If is still close 
to unity, the star should still feel the black hole tidal force but 
will survive the encounter. Hydrodynamical simulations show 
that if the passage is close (0.5 f [5 f \) the star may lose matter 
that can then accrete onto the black hole (Guillochon & Ramirez 
Ruiz 2013). These events should be more frequent than classical 
TDEs (roughly by a factor of ~ 10) and should repeat on the star 
orbital period (if the star is on a bound orbit), increasing their 
observability. On the other side, these events are less energetic, 
because the involved mass captured by the black hole is lower. 

Here, after a brief introduction on total and partial tidal dis¬ 
ruption events (Section 2), we consider the case of a nearby low- 
luminosity AGN, IC 3599, that showed in the past strong X-ray 
activity, possibly resulting from a TDE (Brandt, Pounds & Eink 
1995; Grupe et al. 1995, see Section 3). In Section 4 we describe 
new SwiftpiKf data showing a second, large increase in flux, to¬ 
gether with the analysis of all X-ray data on the source. X-ray 
spectral analysis is presented in Section 5. Based on spectral re¬ 
sults, we converted count rates into a flux and then in luminosity. 
The fit of the overall, 24 yr long, X-ray light curve is described 
in Section 6. Optical data are discussed in Section 7. In Section 
8 we apply partial TDE modelling to the light curve of IC 3599 
trying to constrain the involved star and its orbit. Our conclu¬ 
sions are in Section 9. 
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2. Tidal stripping essentials 

A star orbiting too close to a massive black hole will be torn 
apart by the compact object tidal force. This occurs at the tidal 
radius beyond which stellar self-gravity is not able to counteract 
the black hole tidal field and keep the star together. The tidal 
radius is defined as 

n ~ R, (Mbh/M.)'/ 3 ^ 23 (MiC)-2/3 

where R, and M, are the radius and mass of the star, Mbh the 
mass of the black hole and ~ 3 x 10 ®Mq the mass of the 
central black hole in IC 3599 (Grupe, Komossa & Saxton 2015, 
see also below). The numerical value refers to a IMq and IRq 
star and units are gravitational radii (rg). 

As the star is ripped apart by the tidal forces of the black 
hole, the debris is thrown into high-eccentricity orbits with a 
large range of periods and with an energy range of (Lacy, Townes 
& Hollenbach 1982); 

AC ^ R* 

C 

The distribution of mass as a function of energy is nearly con¬ 
stant (Rees 1988), as also shown by numerical simulations 
(Evans & Kochanek 1989; Lodato, King & Pringle 2009). 

For a parabolic orbit, nearly half of the debris is unbound 
leaving the system at high velocity. The other half will return to 
the black hole at different times depending on the initial eccen¬ 
tricity (i.e. energy). The return time of the first debris at pericen- 
tre is 

2 TT G A/bh 2 tt 

- (2A£)3/2 ~ (GMbh)1/2 (2R.)3/2 “ 

= llRl'^ A 

The return of material at pericentre continues at a rate driven by 
Kepler’s third law as 

3 ^min V^min/ 

where the maximum mass inflow rate is 

Mpeak 1.7R;3/2^2(^IC)-l/2^^yj-l 

‘min 

The peak rate is a function of the stellar structure (Lodato et al. 
2009). After disruption, the stellar debris are launched into very 
eccentric orbits and gradually return to pericentre, where they 
circularise and form an accretion disc at rdrc ~ 2 rp. The follow¬ 
ing fall back of this material onto the black hole is governed by 
viscous times. If the viscous time is short in comparison with 
fmin, then the fall back of matter onto the central object is al¬ 
most instantaneous and fft, ~ fmin- Based on an a— viscosity disc 
prescription we can evaluate the viscous time-scale ty 

^Kep(2 ^p) 

'' nah^ 

where fKep is the Keplerian time at a given radius and h is the disc 
half-height divided by the radius. Being the initial luminosity 
close to (or even somewhat larger than) the Eddington limit, the 
disc is expected to be thick (Ulmer 1999), with h ~ Thus, we 
can estimate the ratio of the viscous time to the time as 

— ~ 1.8 X 

tmin V 10 / 


and a thick disc will form and drain as the material circularises 
down to the last stable orbit. 

The infalling star is usually assumed to be in a parabolic or¬ 
bit and to undergo complete tidal disruption. Encounters may 
still happen with the pericentre slightly outside the tidal ra¬ 
dius (J3 ^ 1), resulting in some spill over of stellar matter, and 
orbits can be bound and highly eccentric, resulting in peri¬ 
odic outbursts. This problem received less attention than the 
classical TDEs but still a number of works exist (Guillochon 
& Ramirez Ruiz 2013; Macleod, Guillochon & Ramirez-Ruiz 
2012; Macleod et al. 2013; Hayasaki, Stone & Loeb 2013; 
Ivanov & Novikov 2001). If the orbit is bounded the physics 
is similar to the parabolic case in terms of the effects of black 
hole tides on the approaching star, especially if the eccentric¬ 
ity is high. On the long term the star orbit may change but the 
typical orbital binding energy is larger than the stellar binding 
energy and the orbital angular momentum is larger than the ro¬ 
tational angular momentum of the star, indicating that the trans¬ 
fer of mass will not substantially alter the orbit (Macleod et al. 
2013). 

In particular, Guillochon & Ramirez Ruiz (2013) investi¬ 
gated the fate of a star undergoing partial stellar stripping. 
They considered the case of stellar polytropes of y = 4/3 
and y = 5/3 and investigated the case for /? within the in¬ 
terval O.b^.O and 0.5-2.5, respectively. They found that to¬ 
tal disruption occurs for a critical impact parameter different 
from unity, being j3c = 1.85 for y = 4/3 and Pc — 0.90 
for y = 5/3, respectively. They also provided fitting formu¬ 
lae for estimating the mass accretion peak rate, the peak time, 
and the total mass transferred for both cases, which are strong 
functions of p, especially for P P Pc (see also the related errata 
corrige at http://astrocrash.net/2013/09/16/errata-of-guillochon- 
ramirez-ruiz-2013/). We used these formulae in deriving stellar 
parameters. 

3. IC 3599 

One of the first putative TDEs occurred in the close active galaxy 
IC 3599 (92 Mpc, at a redshift of z = 0.021; Brandt et al. 1995; 
Grupe et al. 1995). It was discovered in the X-ray band as a 
bright, soft source during the ROSAT all sky survey in Dec. 1990. 
Further ROSAT observations (June 1992 to June 1993) found IC 
3599 in a dimmer (~ 100) and somewhat spectrally harder state 
(see Fig. 1). A Chandra observation in March 2002 found the 
source at a similar level (Vaughan et al. 2004). The very large 
flux decrease and the softness of the X-ray spectrum led sev¬ 
eral authors to suggest that this outburst was a TDE (Brandt et 
al. 1995; Grupe et al. 1995; Komossa & Bade 1999; Vaughan, 
Edelson & Warwick 2004), even if IC 3599 is an active galaxy. 
Indeed active galaxies are predicted to host more TDEs due to 
the perturbing presence of the disc, even if it would be more 
difficult to reveal them because of the higher overall emission 
(Komossa 2012; Karas & Subr 2007). The optical spectrum is 
characterised by strong Ha, HyS and [GUI] lines, showing varia¬ 
tions in response to the X-ray outburst and classifying IC 3599 
as a type 1.5-1.9 Seyfert galaxy (Brandt et al. 1995; Grupe et al. 
1995; Komossa & Bade 1999). The central black hole mass has 
been derived from a relation between the flux at 5100 A and the 
width of the quiescent Hp emission line to be ~ 3 x 10 ®Mq 
(S ani et al. 2010). Grupe et al. (2015) argued that this estimate 
has been obtained using a Broad Line Region scaling relation, 
which is appropriate only in case of an unabsorbed line of sight. 
Using instead a black hole mass to bulge K-band luminosity 
(Marconi & Hunt 2003) or the relation between the [OIII] ve- 
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locity dispersion and the black hole mass (Nelson 2000), they 
found to be in the (2 - 12) x 10 ®Mq range. Here we adopt 
a black hole mass of 3 x 10 ®Mq. For this mass the source was 
well below the Eddington limit during the ROSAT all sky survey 
observation. 

The (low-luminosity) AGN nature of IC3599 may cast 
doubts on the tidal event interpretation, because AGNs show 
flares, which are related to disc activity or to the uncovering of 
a heavily absorbed X-ray source. We note here that the quies¬ 
cent Chandra spectrum is well described by a soft power law 
with minimal intrinsic absorption, weakening considerably the 
absorbed source case. At the same time variations by a factor 
of i 20 - 30 in AGNs are very rarely observed. We study in 
details AGN variability and, based on a sample of highly vari¬ 
able AGNs intensively monitored by Swift, we can assess that 
a variability similar to the one observed in IC 3599 occurs by 
chance at ~ 4.5 cr level (see Appendix A). In particular, the sec¬ 
ond flare should result by chance from known AGN variability 
at i 4.0 cr. Other models apart from AGN variability, have been 
proposed but should not work in the case of IC 3599. A binary 
system made by the central black hole and an orbiting star filling 
its Roche lobe at periastron can be ruled out as the stellar den¬ 
sity required would be too low (Lasota et al. 2011). In addition, 
we confirm that the X-ray and optical transient position is con¬ 
sistent with the centre of IC 3599 (see Appendix B). Periodic 
optical outbursts were observed in the BE Lac object OJ 287 
and explained by accretion instabilities onto a binary black hole 
(Tanaka 2013). However, this mechanism does predict no X- 
rays in quiescence because the innermost part of the disc is de¬ 
pleted. AGN instabilities in slim disk has also been put forward 
to explain this variability (Honma, Matsumoto & Kato 1991). 
However, the estimated mass accretion rate for IC 3599 falls well 
below the instability region and the duration of the flares is much 
longer than theory predictions (Xue et al. 2011, see Appendix). 
All these findings lend support to the idea that IC 3599 under¬ 
went TDEs. Two different TDEs appear however unlikely. The 
occurrence probability of a TDE in a galaxy is of the order of 
10“® galaxy^' yr“^ (e.g. Donley et al. 2002), so having two dif¬ 
ferent events in 25 yr results in a probability of ~ 6 x 10“*. 
Even increasing the rate by a factor of 100 as a result of cen¬ 
tral binary black hole merger (e.g. Perets, Hopman & Alexander 
2007), would result in a still low probability ~ 6 x 10^“*. A bi¬ 
nary disrupted by the central black hole would result into two 
different events but the time delay between them is way too short 
(several days) to account for what we observe (Mendel & Levin 
2015). 

IC 3599 was also detected in the radio band by the Very 
Large Array on June 2012. This radio observation is difficult to 
be accounted for challenging the prediction of the standard blast 
wave model (Bower et al. 2013). 


4. X-ray data preparation 

X-ray data were collected over a time basis of 24 years with 
different satellites. A log of all the observations is shown in Table 
1 . 

4. 1 . ROSAT All Sky Survey data 

RASS data are not straightforward to analyse. They were col¬ 
lected using the PSPC detector in scan mode. Eor this reason we 
decided to stick to the original analysis by Grupe et al. (1995). 


4.2. ROSAT PSPC 

Data were extracted using the task XSELECT (v. 2.4c) from a cir¬ 
cular region centred on source with a 60 arcsec radius. The back¬ 
ground spectra were extracted from an annular region centred 
on IC 3599, free of contaminating sources and support shadows, 
with inner and outer radii 125 and 250 arcsec, respectively. Data 
were retained in the 12-211 channel range, corresponding to an 
energy range of 0.1-2.4 keV. Being IC 3599 always on-axis we 
used the pre-canned response matrix relative to PSPC-B, gain- 
2 period (pspcb_gain2_256.rsp). Given the relatively low num¬ 
ber of photons we binned the data to a minimum of one photon 
per energy channel using the grppha tool and adopted Cash- 
statistics to fit the data. Different observations were grouped to¬ 
gether in order to increase the number of photons as indicated in 
Table 1. 

Cash-statistics does not allow a reliable estimate of the good¬ 
ness of fit. We estimate the goodness of fit by finding the best 
fit with Cash-statistics and then test the data with statistics, 
corrected for the low number of photons per bin through the 
Churazov weighting scheme (Churazov et al. 1996). This allows 
us to determine the goodness of fit using statistics. 

4.3. ROSAT HRI 

Data were first combined into a single event file using XSELECT. 
Source photons were extracted from a circular region of 15 arsec 
radius. Background photons were extracted from a nearby circu¬ 
lar region of 45 arcsec radius. No spectral data can be obtained 
from ROSAT HRI data and the count rate was converted into a 
ROSAT PSPC count rate using PIMMS (see below). 

4.4. Chandra ACIS-S 

Data were reprocessed using the CIAO 4.6 (and CALDB 4.5.9) 
repro task. Spectral data were extracted using the CIAO task 
specextract using a circular 3 arcsec region for the source 
and an annular region for the background with 12 and 17 inner 
and outer radii, respectively. The task provides the user with the 
corresponding response matrix (rmf) and ancillary response file 
(arf). Data were retained in the 0.3-10 keV energy range and 
were binned to a minimum of one photon per energy channel 
using the grppha tool. 

4.5. XMM-Newton Slew survey 

The field of IC 3599 was observed during one XMM-Newton 
slew on 2007 June 19 (MJD 54270). The source was not detected 
with a (2cr) upper limit of 0.8 counts s ' in the EPIC-pn instru¬ 
ment (Grupe et al. 2015). This implies an upper limits on the 
0.3-10 keV unabsorbed flux of 1.4 x 10“*^ erg cm“^ s“* assum¬ 
ing the quiescent spectral model (see below). We do not consider 
this observation any further. 

4.6. Swift XRT 

Data were reprocessed using xrtpipeline (v. 0.12.9) (CALDB 
2014-07-30). Data were extracted from a circular radius of 71 
arcsec and the background from a nearby circular region of 141 
arcsec radius. Data were retained in the 0.3-10 keV energy range 
for spectroscopy and were binned to a minimum of one photon 
per energy channel using the grppha tool. 
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Table 1. IC 3599 X-ray observation log. 


Instrument 

Obs. ID. 

Obs. start 

Duration (ks) 

Counts* 

Spectrum 

ROSAT RASS 

rs931231n00 

1990-12-10 

0.8 

- 

- 

ROSATPSFC 

rp700552a00 

1991-12-15 

3.3 

101 ± 11 

1 

ROSATPSPC 

rp700552a01 

1992-05-31 

1.4 

33 ±6 

2 

ROSATPSPC 

rp701098n00 

1992-06-16 

1.8 

36 ±6 

2 

ROSATPSPC 

rp701099n00 

1992-06-17 

1.9 

33 ±6 

2 

ROSATPSPC 

rp701100n00 

1992-06-18 

2.0 

36 ±6 

2 

ROSATPSPC 

rp701097n00 

1992-06-30 

5.4 

116± 12 

3 

ROSATPSPC 

rp701528n00 

1993-06-17 

3.4 

26 ±5 

3 

ROSAT HRl 

rh702704n00 

1996-06-30 

14.7 

77± 11 

3 

ROSAT ma 

rh702706n00 

1996-06-30 

16.9 

87± 12 

3 

Chandra ACIS-S 

acisf029999N003 

2002-03-07 

10.6 

427 ± 21 

4 

Swift XRT 

00037507001 

2010-02-25 

2.1 

432 ± 23 

5 

Swift XRT 

00037507003 

2010-05-17 

1.2 

100 ± 11 

6 

Swift XRT 

00037507004 

2013-10-30 

4.8 

21 ±6 

7 

Swift XRT 

00037507005 

2013-11-06 

4.9 

12 ±4 

7 

Swift XRT 

00037507006 

2014-03-26 

4.7 

14 ±4 

7 

Swift XRT 

00037569001 

2014-08-08 

1.6 

< 16 

- 

Swift XRT 

00037569001 

2014-11-15 

4.6 

18±5 

7 

Swift XRT 

00037569001 

2014-11-23 

0.4 

< 38 

- 


* The total number of counts was determined using the XIMAGE (v. 4.5.1) task detect (correcting for vignetting and point spread function 
losses). 

Spectra were grouped according to the scheme here reported: spectra with the same number were grouped together. Numbers in italic means that 
we used that given conversion factor to pass from count rates to fluxes even if no spectra analysis were carried out on these data due to the too low 
number of counts. 

Two ROSATPSFC observations on 1992-06-30 and 1992-12-08 were not considered being IC 3599 too far off-axis to be detected. 
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Fig. 1. X-ray spectra fitted with the absorbed disc black body 
plus power law model described in the text. Data were rebinned 
to reach a 5 cr significance or by a factor of 5 for plotting pur¬ 
poses only. In the upper panel black, red and green (filled cir¬ 
cles) data refer to ROSAT spectra (first, second and third observa¬ 
tions, respectively). Blue (open circles) data refer to the Chandra 
spectrum. Light blue, magenta and orange (filled stars) to Swift 
spectra (first, second and third observations, respectively). In the 
lower panel residuals in terms of x are plotted with the same 
colour codings and symbols. 


5. X-ray spectral analysis 

We fitted together three ROSAT PSPC, one Chandra ACIS-S 
and three Swift XRT spectra (see Table 2) with the X-ray spec¬ 
tral fitting package XSPEC (v. 12.8.Ig). All spectral fits were 
minimised using C-statistics and the goodness of the fit was as¬ 
sessed using the Churazov-weighted x^ statistics. It is readily 
apparent that the overall X-ray spectra are soft, as testified by 
previous X-ray data analysis (Brandt et al. 1995; Grape et al. 
1995; Vaughan et al. 2004). An absorbed (using TBABS) power 
law model with all the column densities tied together and the 
same photon index for all the observations returns a;^^^-statistic 
value of 730.9 with 410 degrees of freedom. This power law fit 
provides a null hypothesis probability of 10“^°. The power law 
photon index is very soft L = 3.9 + 0.2 (all errors were deter¬ 
mined for AC = 1.0, i.e. 1 cr errors for one parameter of inter¬ 
est). Leaving free the power law photon index in each spectrum 
improves the fit. We obtained = 493.3 with 404 degrees of 
freedom. The corresponding null hypothesis probability is still 
0.2%. In addition, the photon index of the first ROSAT spec¬ 
trum is extremely high with L = and the column density 

Nh - (4.4 + 0.1) X 10^° cm“^ is much larger than the Galactic 
column density of Nh - 1.2x10^** cm“^ (Kalberla et al. 2005). A 
free column density model with a power law photon index fixed 
instead provides;^^^ = 619.4 with 404 degrees of freedom, corre¬ 
sponding to a null hypothesis probability is 10 ". A soft model 
such as a single black body provides similar results: x^ - 551.7 
with 404 degrees of freedom, with free temperature and radius. 
The null hypothesis probability is 10“^. It is apparent from the 
fit that the black body model fails to account for a high energy 
tail (> 2 keV) at low fluxes. 

We also considered a two component model. Given the large 
number of possibilities we explored two different models, based 
on the physics of possible emission mechanisms. The first model 
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comes from the proposal that this bursting activity comes from 
the partial disclosing of a heavily absorbed AGN. We thus model 
the X-ray spectra with a Galactic absorption plus a partial cover¬ 
ing of a power law component. The second model is physically 
motivated by a tidal disruption event. It consists of an accre¬ 
tion disc spectrum (accounting for the disc emission caused by 
the tidal spoon-feeding) and a power law accounting for the the 
AGN activity (MacLeod et al. 2013). 

Spectral fitting results are shown in Table 2. We explored a 
partial covering factor model (pcfabs within XSPEC) trying to 
see if the outburst observed in IC 3599 can be explained as the 
unveiling of a heavily absorbed source. This model has been en¬ 
visaged to explain the wild erratic variations observed in WPVS 
007 (Grupe et al. 2013). The fit is carried out with a model made 
by a fixed Galactic column density plus a partial covering model 
with the equal intrinsic column densities tied to all the spec¬ 
tra but variable fractions. This composite absorption component 
screens a power law component with the same photon index but 
variable normalisation. We obtained a - 627.4 for 304 de¬ 
grees of freedom, resulting in a null hypothesis probability of 
10-24 (see Table 2). 

In the second model parameters are tied among different 
spectra. The absorbing column density is considered equal for 
all the spectra. The diskbb temperature is expected to vary, 
whereas its normalisation was tied between the spectra (being 
related to the disc inner radius and, likely, corresponding to 
the innermost stable orbit). The power law photon index was 
tied between spectra but its normalisation was free to vary. The 
overall fit provides a x^ = 315.9 for 303 degrees of freedom. 
The null hypothesis probability is 29%. Residuals are well dis¬ 
tributed all over the entire energy band (see Fig. 1 and Table 
2). The derived column density is somewhat in excess of the 
Galactic value. From the normalisation of the disc black body 
model one can estimate the inner disc radius (modulo the disc 
inclination, i, and a colour factor uncertain by a factor of fc ~ 2). 
Considering a maximally rotating black hole, we can estimate it 
mass as ~ 8 X 10^ x (cos 112)^1^ x (/c/2)Mo. This estimate, even 
if approximate, is barely consistent with the value derived from 
optical studies, indicating that the central black hole in IC 3599 
is not particularly massive. 


Table 2. IC 3599 X-ray spectral fits. 



diskbb-l-pow 

pcfrpow) 

Nh (10^") cm-^ 

PLT 

Diskbb norm. 

O Q+U.2 

1900+^*^ 

4 4 + 1-0 

3648:? 

Dataset 

r(eV) 

Cov. fraction 

ROSl 

ROS2 

ROS3 

Chal 

Swil 

83« 

75/4 

77+4 

ic\ 

142+* 

116/2 

034+0177 

^•^^- 0.07 

^•^^- 0.02 

0.95/°;°! 

0.80/|i;i 

<0.10 

Swi2 

< 0.17 

Swi3 

51^^ 

-’'-16 

0 984° °' 

(dof) 

315.9 (303) 

627.4 (304) 

nhp 

0.29 

1 

o 


Errors were determined with AC = 1. 


6. Light curve fitting 

We fit together the luminosity light curve derived from all Swift 
X-ray data and archival ROSAT and Chandra data. To do this 
we computed the conversion factors from count rates to unab¬ 
sorbed fluxes in the 0.01-10 keV energy band for each spectrum 
based on the absorbed disc plus power law spectral model. These 
fluxes were then converted into luminosities adopting a distance 
of 92 Mpc for IC 3599 and used in Fig. 2. Spectra that were ob¬ 
tained as the sum of different observations were then split into 
single observations and the same conversion factor was applied 
to all of them. The flux of last ROSAT observation was es¬ 
timated based on the spectrum of the third ROSAT PSPC obser¬ 
vation (i.e. the closest in time). The flux of the RASS point was 
extrapolated to the 0.01-10 keV energy band based on the best 
fit spectrum (Grupe et al. 1995). For the ROSAT HRI count rate 
we adopted a different approach. We fitted the last ROSAT PSPC 
spectrum with absorbed single component models (power law, 
black body and bremsstrahlung). Based on these models we used 
PIMMS (v. 4.7b) to extrapolate the ROSAT HRI rates to equiv¬ 
alent ROSAT PSPC rates. An average rate was computed and 
the standard deviation among the different models was added in 
quadrature. Then we converted these rates into fluxes as above. 

The overall Swift light curve shows many similarities with 
the previous ROSAT light curve (Brandt et al. 1995; Grupe et 
al. 1995; Komossa & Bade 1999). Both bolometric light curves 
can be fitted with a power law, pointing to (at least) two 
different tidal disruption events. The similarities in the decay 
of the events suggest that we are observing a recurrent phe¬ 
nomenon rather than the random occurrence of different accre¬ 
tion episodes. The non-disruptive passage of the same star pro¬ 
vides a more comprehensive explanation. 

We started by fitting the luminosity light curve with a differ¬ 
ent number of outbursts. Values of the x^ reported in Table 
3. It is readily apparent that the fit with 3 outbursts and a power 
law shape is superior to all the others. The improvement of a free 
power law index with respect to the fixed value -5/3 has been 
evaluated by means of an F-test. The probability of a random 
occurrence is 12%. 

In addition to the luminosity light curve, we folded the tem¬ 
peratures derived from spectral analysis along with the sug¬ 
gested orbital period. Disc temperature should evolve as T oc 
j-5/12^ being M oc and T oc mF4 (from accretion disc the¬ 
ory, Lodato & Rossi 2010). Fitting the temperature evolution for 
three outburst peaks we obtain a reduced x^^^ - 0.6 for 6 de¬ 
grees of freedom and a null hypothesis probability of 75% (up¬ 
per panel in Fig. 3). If we fit instead the temperature evolution 
for two outbursts only with the same model we derive a reduced 
X^^^ = 32.0 for 6 degrees of freedom and a null hypothesis proba¬ 
bility of 10-2* (orange dashed line in Fig. 3). If we add a constant 
to the power law model, mimicking the presence of an underly¬ 
ing quiescent accretion disc we obtain= 6.6 for 5 degrees of 
freedom and a null hypothesis probability of 4 x 10-® (continu¬ 
ous line in the lower panel of Fig. 3). The constant quiescent disc 
temperature is 42 + 6 eV. Note that such a constant temperature 
baseline might be present also in the first fit but it is not required 
by the data. We note that the presence of three peaks depends 
on the Chandra data only. If these data are strongly affected tim¬ 
ing and spectral variability unrelated to the TDF, a longer period 
(twice) should be envisaged. 

Based on light curve and disc temperatures evolution we 
found strong indication for multiple (equal) tidal stripping 
events. In particular, the timing and spectral data points toward 
three (rather than two) outbursts taking place during the 1990- 
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Time (MJD) 


Fig. 2. Long term X-ray luminosity light curve of IC 3599. Error bars are at 1 cr confidence level. ROSAT, Chandra and Swift count 
rates were converted into 0.01-10 keV unabsorbed luminosities by means of spectral fits assuming a source distance of 92 Mpc. The 
overall X-ray light curve has been ht with a (f - fo) function repeating over a Pq time. Free parameters are the starting date, the 
power law normalisation, the repetition time and a constant, setting the basic emission level of IC 3599. The quiescent level (AGN 
activity) has been evaluated worsening the ht until a 10% null hypothesis probability is attained. This is nicely consistent with the 
mean value and the observed variability in the power law component used in the X-ray spectra hts. 


Table 3. X-ray luminosity light curve htting. 



(dof) 

nhp 

r“ 

(dof) 

r“ 

nhp 

a 

Exp. 

^2 (dof) 

Exp. 

nhp 

e- fold 
(d) 

2 outburst 

37.3 (12) 

2 X 10-^ 

25.3(11) 

8 X 10-3 

-3.8+ 1.2 

33.0(11) 

5 X 10-^ 

108 + 9 

3 outburst 

9.3 (12) 

7 X 10^* 

7.1(11) 

8 X 10-‘ 

-2.7+ 1.1 

28.4(11) 

3 X 10-3 

112 + 6 

4 outburst 

29.0 (12) 

4 X 10-3 

24.3(11) 

1 X 10-2 

-3.6+ 1.6 

33.0(11) 

5 X 10-^ 

108 + 6 


2014 time span, with a recurrence time of 9.5 yr (see Figs. 2 and 
3). Assuming an accretion efficiency of 77 = 0.1, from the out¬ 
burst light curve we can estimate the peak mass accretion rate to 
be Mpeak ~ O.OIMq yr“* and the accreted mass per episode to 
be AM ~ 2.5 x IO^^Mq. Clearly these are lower limits, missing 
the early stages of all the outbursts. In addition, with these new 
ephemerides, radio observations took place only 2.6 yr after an 
outburst episode and can be accounted for by emission internal 
to the jet (van Velzen, Kording & Falcke 2011). 


7. Optical data 

IC 3599 fell within the field of view of the Catalina sky survey 
and was monitored in the optical in the MJD 53470-56463 inter¬ 
val (Drake et al. 2009). The Catalina light curve shows a broad 
peak (see Fig.|^. The light curve can be fit with a Gaussian cen¬ 
tred on MJD 55151 + 11 (Icr) and with width of 200 ± 18 d. The 
peak of the optical emission occurs after the estimated peak of 
the X-ray emission by ~ 140 d. However, due to the large width 
of the optical flare, the optical emission starts ~ 280 d before the 
high energy peak (assuming as the start of the optical flare when 
the optical flux rises by 10% over the constant value). 
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Time since outburst peak (d) 


Fig. 3. Upper panel, disc temperatures as derived from the spec¬ 
tral fits folded on a three passages light curve. ErTor bars are at 
1 cr confidence level. Red points (ROSAT) refer to the first pas¬ 
sage, the green point (Chandra) to the second passage and blue 
points (S wift) to the third one. Disc temperature evolution is ht- 
ted with a hxed T oc power law (deriving from M oc 
and T oc typical of an accretion disc and therefore be¬ 

ing a direct test for the disc cooling, Lodato & Rossi 2010). 
Lower panel, as above but in the case of two passages. Red dots 
{ROSAT and Chandra) refer to the hrst passage and green dots 
(S wift) to the second one. The ht with a T oc power law is 
shown as a dashed orange line. 

This gives us the possibility to study in detail the outburst 
start (and prepare a follow-up strategy for the next outburst). The 
rise time in a TDE is governed by the circularisation time (i.e. 
the time to form the accretion disc) and the viscous time (i.e. the 
time needed to transfer matter from the outer disc edge to the 
central compact source). Several papers have appeared recently 
on this subject (Bonnerot et al. 2015; Guillochon & Ramirez- 
Ruiz 2015; Piran et al. 2015). Irrespective of which is the domi¬ 
nant mechanism, we can safely assume that the viscous time of 
the disc, C, is shorter than the observed rise time. Given the opti¬ 
cal light curve the time it takes to rise from 10% of the quiescent 
flux to the peak flux is ~ 420 d. Eitting the optical light curve 
with a symmetric exponential function, we derive an e-folding 
rise/decay time of 178 + 14 d. Based on Guillochon & Ramirez 
Ruiz (2013), the viscous time can be expressed as: 

C ^ 12A 


8. Constraining the orbit and star characteristics 

Partial disruption of stars has been investigated by means of hy- 
drodynamical simulations (Guillochon & Ramirez Ruiz 2013; 
Macleod et al. 2013; Macleod et al. 2012). Depending on the 
star structure (modelled as a polytrope with index y - 4/3 or 
5/3), it has been shown that mass can be extracted from the star 
for impact parameter (3 as low as 0.5-0.6 whereas complete de¬ 
struction occurs for (3 > 0.90 for 7 = 5/3 and for (3 > 1.85 for 
y - 4/3, respectively (Guillochon & Ramirez Ruiz 2013). In 
order to constrain the encounter and star characteristics we com¬ 
pare the mass accretion at peak rate Mpeak and the total amount 
of matter 2 AM lost by the star during each passage with the cor¬ 
responding quantities estimated through hydrodynamical simu¬ 



Fig. 4. Optical (flux) light curve of IC 3599 from the Catilina sky 
survey in arbitrary units (red dots), together with the SwiftpCKY 
light curve (open squares) converted in flux using the spectral fit 
model, in arbitrary units. The Gaussian fit refers to the optical 
data. The vertical dashed line marks the time (MJD) of the op¬ 
tical peak, the dotted line the peak of the X-ray flux according 
to the multi-peak modelling. The continuous line marks the time 
at which the optical light curve rises by 10% from the quiescent 
level. 


lations (Guillochon & Ramirez Ruiz 2013). These quantities de¬ 
pend strongly on the impact parameter j3 and, to a lesser extent, 
on the star mass and radius. 

The peak mass accretion rate has been estimated based on 
the RASS data, converting the peak luminosity into a peak mass 
accretion rate by assuming a 10% conversion efficiency. In order 
to bracket uncertainties the allowed parameter space is computed 
accepting values of the peak accretion rate in the interval 1-3 
of observed value. The other parameter considered is the total 
amount of mass accreted during one flare episode. This has been 
estimated by integrating the luminosity light curve over one or¬ 
bital period. Also for this parameter we searched solutions in the 
interval 1-3 of the observed total mass accreted. Finally, we re¬ 
quire that the viscous time of the disc, ty, is shorter than 178 + 42 
d (3 cr). A mass range of 0.2 - IOOMq and a (unrelated) radius 
range of 0.01 - 10"^Ro were blindly searched. We consider two 
different star models based on a polytropic index of y = 4/3 
(over the interval (3 - 0.5 - 4.0) and y = 5/3 (/? = 0.45 - 2.5). 
Results are shown in Fig. 4 for the two indexes separately. 

In case of y = 4/3 we And solutions in the range (3 - 
0.58 - 0.72. The range of allowed values of (3 has been obtained 
without requiring a priori a partial disruption event, searching 
for solution in the (3 - 0.5 - 4.0 interval. This results therefore 
strongly (and independently) supports that the flaring events in 
IC 3599 are related to a partial tidal disruption (y6 < = 1 -85) of 

an orbiting star. The corresponding mass range is for M, > 4Mo. 
We blindly explored the mass-radius plane, but stars do not All 
this plane homogeneously. To test the consistency of our find¬ 
ings with stellar models, we investigated the mass-radius relation 
by means of the Single Stellar Evolution (SSE) code (Hurley, 
Pols & Tout 2000), assuming solar metallicity (we verified that 
a change in the metallicity does not change our results sensi¬ 
bly). With this additional constraint we have that the range of 
allowed masses reduces to M» = 15 - 45 Mq for stars on the 
main sequence and to M* =4-15 Mq for slightly evolved stars. 
The corresponding radii in the selected region are R, ~ 5 - 9Rq 
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for main sequence stars and R, ~ 4 -1 Rq for slightly evolved 
stars (see Fig. 4 upper panel). The allowed eccentricities range 
in e ~ 0.995 - 0.997 and the pericentre distance rp ~ 84 - 93 
(where is the gravitational radius of the black 

hole in 1C 3599, G the gravitational constant and c the speed of 
light). The eccentricity is high but not unlikely given the pre¬ 
dicted distribution of eccentricities around the central black hole 
in our Galaxy (Gillessen et al. 2009). This eccentricity is larger 
than the critical eccentricity for bound orbits below which all 
stellar debris remains bound and feed the black hole on a longer 
time scale (Hayasaki et al. 2013). 

In the case of a poly tropic star with index y = 5/3, we 
searched solutions for fi in the 0.45-2.5 range (see Fig. 4 lower 
panel). We do hnd solutions in a narrow range ofjS = 0.49-0.55. 
Again the allowed range of jS is well below the critical value for 
a complete tidal disruption {Pc - 0.9). The allowed mass range 
is M* 1.5Mo. However polytropic stars with y = 5/3, can sat¬ 
isfactorily describe Sun-like stars (or smaller) and we do not hnd 
therefore acceptable solutions in this case. 

9. Conclusions 

We report here on the hrst candidate of periodic, partial tidal dis¬ 
ruption events. These events were observed as large (> 100) hux 
increases in the 24-yr long X-ray light curve of the close, low- 
luminosity active nucleus IC 3599 (see Fig. 2). We discuss sev¬ 
eral possibilities to explain these large hares and likely conclude 
that they cannot be ascribed to AGN variability; single hares are 
unlikely, at i 4.0 cr level, to come from known, observed ex¬ 
tremely variable AGNs, and the haring instability region sug¬ 
gested by some accretion disc models (e.g. Honma et al. 1991) 
lies above the mean accretion rate observed in IC 3599 by a fac¬ 
tor of i 15. By modelling the light curve and the disc tempera¬ 
ture evolution, as derived by spectral hts, we hnd out that three 
Hares are statistically preferable over two, resulting in a ~ 9.5 
yr periodicity. This cannot be appreciated from the light curve 
where the second Hare is largely missed due to sparse observa¬ 
tions, but results from the disc temperature evolution: if the stan¬ 
dard disc model can apply to the observations of IC 3599, then 
the 3-outbursts model is preferable over the 2-outburst model at 
~ 3.8 cr. Based on the disc temperature evolution the same is true 
at ~ 4.3 cr conhdence level. A tidal disruption model applied to 
IC 3599 provides a good match to the data and, independently, 
solutions for partial disruption events only. 

In addition, we note that the average quiescent luminosity of 
the AGN is just a fraction of average Hare luminosity (~ 20%), 
leaving open the possibility that the overall AGN activity of IC 
3599 is entirely spoon-fed by the orbiting star (MacLeod et al. 
2013). Given the relatively large mass stripped every passage, 
the phenomenon is short-lived and based on our estimates it can 
last .iS 10'^ yr. The next passage should occur in 2019, giving us 
the ability to plan a detailed monitoring campaign to explore the 
characteristics of the orbiting star. The very short orbital period 
gave us the possibility to observe more than one passage. Other 
tidal disruption events presently known might be similar to this 
case but just with a longer orbital period. 
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Appendix A: AGN variability 

Given the low-luminosity AGN nature of IC 3599 we investi¬ 
gated in details here if the observed flares can come from this 
activity. We approached the problem either studying standard 
AGN variability and by taking an unbiased sample of the most 
variable AGNs as observed by the Swift satellite. 

A. 1. Intrinsic AGN variability 

We simulated the long-term light curve of an AGN (Vaughan 
et al. 2003). We assumed a standard broken power law for the 
power spectral density (PSD) of a long term (years) monitor¬ 
ing light curve. We considered an index -1 at short frequencies 
and -2 at high frequencies, with a break at 10^ Hz scaled to the 


mass and luminosity of the black hole in IC 3599 (McHardy et 
al. 2006). We checked that even heavily changing these num¬ 
bers the final result does not depend on them. We run a sim¬ 
ulation generating a number of light curves based on the PSD 
above, including no root-mean-square (rms) variability and no 
background. The long-term light curve changes in all cases are 
less than a factor of a few with respect to the mean starting value 
(on timescales of days). We then turned on the rms variability 
and this value is the main driver for the variability. Even assum¬ 
ing a 100% rms variability, we obtained a maximum increase by 
a factor of ~ 6 - 7 in the count rates (no background included) 
from the mean value on timescales longer than days. From this 
analysis we conclude that normal AGN variability is not able to 
produce the strong variability we observe. 

A.2. AGN with large intrinsic variations 

We started from SwiftfXRT light curves spanning a 10 yr basis. 
We first considered the AGN sample monitored by BAT, consist¬ 
ing of bright nearby AGNs with more than 5 SwiftpCRT observa¬ 
tions. In this sample there are just 3 AGNs showing a count rate 
change by more than a factor of 20. These are Mkn 335 (A = 42), 
NGC 4395 (A = 141) and CGCG 229-015 (A = 56, even if 
in this case the change is a dip rather than a flare). One very 
efficient method to And highly variable AGN consists in com¬ 
paring soft X-ray observations at different epochs. Grupe et al. 
(2001) showed that 4 out of 113 bright RASS objects observed 
also in ROSAT pointings showed dramatic variations. Grupe et 
al. (2012) added other 4 objects to this list by means of subse¬ 
quent XMM-Newton and Chandra observations. All but two are 
Narrow Line Seyfert 1 galaxies. The first is RX J1624.9-I-7554, 
a non-active galaxy for which a TDE has been invoked (Grupe 
et al. 1999), the second one is IC 3599. In this Grupe-Komossa 
sample of highly variable AGNs there are then 8 objects but only 
3 of them show variations larger than a factor of 30. These are IH 
0707^95 (A = 114, dip-like variation), Mkn 335 (as above) and 
WPVS 007 (A = 36). WPVS 007 is however too faint for Swift 
and the light curve is unusable because it has too few points 
(Grupe et al. 2013). We have then 4 objects with good-enough 
sampled light curves, showing flare-like features (Fig. A.l). 

One conservative test is to extract randomly from these 
curves 7 points (the number of points comprising the first TDE 
flare) and 6 points (the number of points comprising the third 
TDE flare) separately, and fit them with a tidal decay template 
(f - -H c to mimic one TDE at a time. We simulated 500,000 
light curve realisations for each source by randomly sampling 
the observed light curves, extracting 7 or 6 points, respectively. 
We then excluded simulated light curves with a rate variation 
less than a factor of 30 (in order not to fit constants, note that 
even if the flux variation in the Swift/XRT data is ~ 30 the count 
rate variation is ~ 100 due to spectral variability, therefore our 
rate change provides a conservative estimate) and fit the remain¬ 
ing light curves with the above model. We adopted two different 
approaches. In the first one we counted how many simulations 
have a null hypothesis probability larger than 5% and, based on 
this number, evaluated the probability of obtaining by chance a 
TDE-like light curve. In the second approach we took among all 
these light curves the best one in terms of and derived the 
probability of this with respect to the number of degrees of 
freedom (5 and 4 in our case, respectively). We do this because 
even the best selectedis, for some sources, not good, at vari¬ 
ance with the fit of the IC 3599 flare. This x^ probability is then 
weighted for the number of trials (500,000). With these num¬ 
bers we assessed the probability of randomly extract a TDE-like 
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Fig. A.l. SwiftPiKT light curves of the 4 selected highly variable 
AGN. 


event from these light curves with a count rate increase by a fac¬ 


tor of > 30 (the Swift rate variation is ~ 100). We obtained in 
this way conservative limits that are shown in Table A. 1. 

A different check can be made on the spectral properties of 
IC 3599. In Fig. A.2 we show the spectral evolution of the hve 
sources in our sample. This is quantified in terms of hardness 
ratio (computed as HR - (H - S)I(H -i- 5), being H and S the 
counts in the 2-10 keV and 0.5-2 keV energy bands, respec¬ 
tively, so that HR can vary among -1 and 1) versus cont rates 
normalised to the minimum observed rate. In this plane IC 3599 
stands alone, reaching complete softness {HR - -1) at maxi¬ 
mum and showing a clear spectral evolution. IH 0707-495 is 
similarly soft but it does not show a marked spectral variability 
as the count rate changes. Mkn 355 and NGC 4395 show spectral 
marked changes but never appear as very soft sources. 

We should safely conclude that the flares observed in IC 
3599 do not likely come from known (observed) AGN variabil¬ 
ity. 
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Fig.A.2. Hardness ratio as a function of the normalised count 
rates for the variable AGNs described above. IH 0707^95 data 
are shown as light blue open stars, CGCG 229-015 as black 
open squares, NGC 4395 as blue open circles, Mkn 335 as green 
open triangles, and IC 3599 as filled red circles. Count rates are 
normalised to the minimum rate observed. The hardness ratio is 
computed as the ratio as (H - S)/{H + S) being H the number 
of counts in the 2-10 keV energy band and S those in the 0.5-2 
keV band. 


A.3. Disk instability models 

The standard geometrically thin, optically thick accretion disc 
model has thermally unstable regions if the accretion rate is 
larger than a critical value. A new branch of equilibrium so¬ 
lutions should exist for rate across the Eddington luminosity, 
named slim discs (Abramowicz et al. 1988). Some works have 
shown that a limit cycle might take place, resulting in a flare- 
like behaviour (Honma et al. 1991). Magneto-Hydrodynamical 
simulations showed, on the contrary, that such disks could be 
thermally stable (Hirose, Krolik & Blaes 2009). Observationally, 
Galactic transient X-ray binaries in outburst do not show any 
evidence for such an instability, except perhaps GRS 1915-1-105 
(Xue et al. 2011). 

The mean mass accretion rate of IC 3599 is ~ 0.005 in 
Eddington units (for an efhciency of 10%, where the peak ac¬ 
cretion rate reaches ~ 0.1 the Eddington rate), far away from 
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Table A.l. Simulated light curves from Swift/XRT variable AGN searching for random TDE-like events. 


Source 

TDEl probab. 

(cr) 

TDE3 probab. 
(cr) 

Num. Obs. 

Variation 

Type 

IH 0707^95 

4.11(4.95) 

4.06(4.84) 

101 

A = 114 

dip 

Mkn 335 

> 4.75(7.93) 

> 4.75(7.89) 

179 

A = 42 

flare 

NGC 4395 

> 4.75(5.74) 

4.06(4.94) 

190 

A = 141 

flare 

CGCG 229-015 

>4.75(> 8.31) 

>4.75(> 8.31) 

44 

A = 56 

dip 


Probability of the first and third IC 3599 flares to be obtained by sampling the Swi/f/XRT light curves of highly variable AGNs. The first probability 
is obtained by requiring that the simulated light curve has a null hypothesis probability larger than 5%. The second probability (in parenthesis) is 
obtained by weighing the best fit light curve probability with the number of simulations carried out. 


border of the unstable region, starting at a critical mass accre¬ 
tion rate > 0.1 in Eddington units (Honma et al. 1991). Elare 
models do not make a clear prediction of the shape of the flare 
but make key prediction that the flare width is a very short frac¬ 
tion of the times between two different flares (Xue et al. 2011). 
The ratio of the EWHM of the flare with respect to the duration 
time between two consecutive flares must be i 3.6%. We mod¬ 
elled the two main peaks of IC 5399 with a Gaussian, limiting 
the maximum luminosity to a few times the Eddington value, 
and found a EWHM~ 500 d (leaving free the maximum lumi¬ 
nosity provides a better ht and an even larger EWHM~ 1200 d). 
Taking the times among the two main flares, we obtain a ratio 
of ~ 7% (~ 17% for a free Gaussian shape). This is more than 
double what expected based on models. 

We also note that two possible very recent flare-like events in 
AGNs, XMM SLl J061927.1-655311 (Saxton et al. 2014) and 
NGC 2617 (Shappee et al. 2014), involve smaller flux variations 
( C 20) and much more complicated X-ray light curves. 

Appendix B: UVOT analysis 

B.1. Image analysis 

Together with XRT exposure, the Swift satellite took contempo¬ 
raneous images of IC 3599 with the UVOT instrument. During 
all the six Swift observations (but the fourth), UVOT observed IC 
3599 with ah its 6 Alters (in the fourth only with the uvm2 Alter). 
Sky images were considered and fluxes were obtained with the 
uvotsource task using the latest calibrations (Poole et al. 2008; 
Breeveld et al. 2011). Eor UV images a circular extraction region 
of 8 arcsec radius has been considered and for optical image a 
5 arcsec radius, respectively. A background region close to IC 
3599 free of sources and of 20 arcsec radius has been selected. 
Results are shown in Table B.l. A clear decrease in flux is ap¬ 
parent, especially at UV wavelengths with a decrease by a factor 
of ~ 2 (uvm2). 

B.2. Astrometry 

In order to determine precisely the position of the transient 
source observed in IC 3599, we carried out digital image sub¬ 
traction between the flrst and the last UV frames (uvm2) ob¬ 
tained with the Swift/UVOT. In the flrst image, obtained on 2010 
Eebruary 25, the transient source is at its maximum while in the 
last epoch, obtained on 2014 March 26, it has faded below the 
host galaxy level (Table B.l). Both images were obtained with 
the uvm2 Alter and have exposure times of 0.2 ks and 1.2 ks, 
respectively. Before subtraction, the two images were aligned. 
We accounted for the difference in the exposure times by mul¬ 
tiplying the image obtained on Eeb 2010 by a factor of six. As 
can be seen in the third panel of Eig. B.l, the result of the im- 


Table B.l. IC 3599 UVOT data analysis. 


Time 

Flux 

(MJD) 

10“'^ erg A“' 

V band 

(5468 A) 

55252.3 

2.23 ±0.11 

55333.0 

2.00 ± 0.09 

56594.9 

1.70 ±0.10 

56741.9 

1.61 ±0.08 

b band 

(4392 A) 

55252.3 

2.20 ± 0.09 

55333.0 

2.03 ± 0.07 

56594.9 

1.63 ±0.08 

56741.9 

1.58 ±0.06 

u band 

(3465 A) 

55252.3 

1.97 ±0.08 

55333.0 

1.79 ±0.07 

56594.9 

0.99 ± 0.06 

56741.9 

1.03 ±0.05 

uvwl band 

(2600 A) 

55252.3 

2.43 ±0.10 

55333.0 

2.20 ± 0.08 

56594.9 

1.09 ±0.06 

56741.9 

1.19 ±0.06 

uvm2 band 

(2246 A) 

55252.3 

2.59 ±0.12 

55333.0 

2.55 ±0.11 

56594.9 

1.16 ±0.08 

56601.9 

1.22 ±0.04 

56741.9 

1.23 ±0.04 

uvw2 band 

(1928 A) 

55252.3 

2.90 ± 0.09 

55333.0 

2.71 ±0.08 

56594.9 

1.43 ±0.06 

56741.9 

1.33 ±0.05 


Errors are at 1 cr confidence level. 

age subtraction is good, showing a clear residual at the position 
of IC 3599. The position of this residual is RA, Dec (J2000): 
12:37:41.18, -t26:42:27.0 (with an uncertainty of 0.3 arcsec, cal¬ 
ibrated against the USNOBl.O catalogue). At the distance of IC 
3599 1 arcsec corresponds to 0.4 kpc. The residual is located at 
0.8 + 0.4 arcsec (1 cr confidence level) with respect to the posi¬ 
tion of the IC 3599 galaxy centre and is fully consistent with the 
position of the candidate TDE radio counterpart (Bower et al. 
2013). The consistency of the position of the UV transient with 
the IC 3599 galaxy centre (at the 2 cr level) is in agreement with 
a TDE origin. 
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Fig. B.l. Swift/UVOT images of IC 3599, obtained with the uvm2 filter on 2010 Feb 25 (left panel) and 2014 March 26 (central 
panel). The images are 5' x 5', North is up, East is left. The right panel shows the result of digital image subtraction between these 
two images. A clear residual is visible at a position consistent with the IC 3599 galaxy centre. 
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